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Abstract 

We introduce an anisotropic two-dimensional lattice gas model of metal terminated 
II-VI(OOl) semiconductor surfaces. Important properties of this class of materials are 
represented by effective NN and NNN interactions, which result in the competition of 
two vacancy structures on the surface. We demonstrate that the experimentally observed 
c(2 X 2)-(2 X 1) transition of the CdTe(OOl) surface can be understood as a phase transition 
in thermal equilibrium. The model is studied by means of transfer-matrix and Monte 
Carlo techniques. The analysis shows that the small energy difference of the competing 
reconstructions determines to a large extent the nature of the different phases. Possible 
implications for further experimental research are discussed. 



Two-dimensional lattice gases have served as models of atoms adsorbed to a singular crystal 
surface, or the terminating layer of such a surface itself, respectively. The interplay of at- 
tractive and repulsive short range interactions can result in highly non-trivial features, see 
e.g. [Q, ^, ^, ^, |5| and references therein. For instance, square lattice systems with infinite 
NN-repulsion (hard squares) and NNN-attraction display tricritical behavior. At low tem- 
peratures a dense, c(2 x 2) ordered phase coexists with a disordered phase of low coverage. 
Here we will investigate a particular model with highly anisotropic attractive and repulsive 
interactions, which result in a c(2 x 2) groundstate, as well. However, this ordering competes 
with a (2 X 1) structure which can prevail locally in the disordered regime. 

The model parameters are chosen as to represent certain properties of metal terminated II- 
VI(OOl) semiconductor surfaces. This class of materials has attracted considerable attention 
due to their potential technological relevance in the development of optoelectronic devices, 
for a recent overview see Frequently, (001) surfaces serve as substrates for the growth of 
II- VI crystals by means of Molecular Beam Epitaxy or Atomic Layer Epitaxy, for instance. 
Surface reconstructions play an important role in this context and have been the target of 
experimental studies ^, In contrast to most III-V materials, II-VI(OOl) surfaces exhibit 
a fairly small number of possible reconstructions, which are less complex than their III-V 
counterparts, in general. 
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In the following we will mainly address the CdTe(OOl) surface, see for a detailed dis- 
cussion. Apparently, only Cd-terminated (001) surfaces are observed in vacuum 0,0]. The 
underlying, complete Te half-layer provides potential Cd-sites which form a simple square 



lattice. Electron counting rules and similar considerations |13] show that the simultane- 
ous occupation of NN-sites in the [1 1 0] -direction (termed the y-direction in the following) 
is excluded in the terminating Cd-layer, whereas NN-neighbors along the [110]-direction 
(or X-axis, for short) are possible. Therefore, unless excess Cd is deposited, the surface is 
characterized by a vacancy structure with a maximum Cd-coverage of ^ = 1/2. 

Figure |l| (a) illustrates the structure of the two relevant configurations which satisfy this 
constraint at 9 = 1/2. The c(2 x 2) reconstruction is characterized by a staggered {checkered) 
occupation of the square lattice sites. In the (2 x 1) structure, Cd-atoms arrange in rows 
along the x-direction which alternate with rows of vacancies. In principle, the configurations 
can be transformed into one another by shifting every other column of Cd-atoms by one 
lattice site. 

Density functional (DF) calculations have shown that the surface energies of the two 
competing structures at 9 = 1/2 and T = differ only by a small amount AE, with the 
c(2 X 2) reconstruction having the slightly lower energy. Qualitatively this preference can be 
unterstood in terms of electron Coloumb interactions, as the distances of neighboring metal 
atoms are smaller in the (2 x 1) arrangement |14]. For ZnSe, a value of AE ~ O.OSeF per 
potential Zn-site is given in [H, m, |l|]. According to 0, the energy difference is even 
smaller {AE ^ 0.016eF) for the CdTe(OOl) surface. 

This factor should play a crucial role in a phase transition which has been studied for CdTe 
1^, |8|, 1^: in vacuum at temperatures below a critical value of about Tc = 270°C ± 10°C, the 
surface displays a mixed c(2 x 2)-(2 x 1) structure with a clear prevalence of the checkered 
configuration close to (but below) Tc- Above Tc, the (2 x 1) arrangement of Cd-atoms 
dominates the surface. The observed coverage is in the vicinity of ~ 0.4 in both regimes 
j^]. The situation is complicated by the fact that the material begins to sublimate at about 
the same temperature Tc- However, it has been argued that sublimation through step flow 
would not hinder the surface to achieve an effective equilibrium configuration on terraces j^. 

The aim of our theoretical investigation is to clarify, whether the nature of the above 
discussed transition can be explained within a thermodynamic equilibrium framework at all, 
or if non-equilibrium effects should play a crucial role. 

The modeling of reconstructions which are characterized by displacement of atoms from 
their regular lattice positons, usually requires continuous two- or three-dimensional degrees of 
freedom. A prominent example is the description of W(IOO) surfaces by XY-models, see e.g. 
[|l8|| and references therein. Here, however, reconstruction occurs via the rearrangement of 
atoms in vacancy structures and a description in terms of occupation variables is appropriate. 

We present here a lattice gas model which takes into account important features of the 
above discussed II-VI(OOl) surfaces. We will loosely speak of Cd-atoms in the following, 
without claiming to reproduce particular properties of CdTe faithfully. In fact, the basic 
structure of the model would be the same for other II-VI(OOl) surfaces. In our simplifying 
picture we consider only the terminating Cd-layer, represented by a square lattice of sites 
{x,y) which can be either occupied {rix^y = 1) or empty {rix^y = 0). The influence of the 
underlying crystal structure is accounted for by effective pairwise interactions of atoms. In 
the 7/-direction, an infinite repulsion excludes the simultaneous occupation of NN-sites, i.e. 
n'x,y = 1 always implies n^^y^ = 0. In the x-direction, an attractive interaction favors the 
occupation of NN-pairs, the strength of which is denoted by J^; < 0. A competing attractive 
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interaction of diagonal neighbors (NNN) < tends to stabilize the c(2 x 2) arrangement 
of atoms. The total energy of the system is given by 

H = "^.y ( [f^X+l,iAfl + "-X+l,S/-l] + Jxn,j;+l^y - flj , (1) 

where the smn is over all lattice sites and the (effective) chemical potential /z controls the 
mean coverage 6 = { n^ y ) < 1/2. Without loss of generality we can choose = —1 and thus 
fix the energy scale. Then controls the energy difference AE (in units of \Jd\) between a 
perfectly ordered c(2 x 2) and a perfect (2 x 1) arrangement at 6 = 1/2: AE = \2 + \/2 
(per lattice site). The groundstate of the system is a c(2 x 2) ordered configuration with 
9 = 1/2, whenever Jx > —2 (and /u > —2). 

The free energy of the system is obtained from the partition function Z = J2{nx y} 
where the temperature T = 1//3 is also measured in units of | | = 1. The sum is restricted 
to configurations { Ux^y } which obey the NN-exclusion in y-direction. We have applied 



standard transfer matrix (TM) techniques [19| to evaluate the logarithm of Zl, the partition 



sum of a system with M = L x N lattice sites in the limit N —)■ oo. Strips of width L with 
periodic boundary conditions were aligned with the x-axis. Hence, only even L allow for the 
perfect c(2 x 2) ordering of the groundstate. Note that the TM is of dimension 2^ x 2^, but 
with a much smaller number 3^ of non-zero elements due to the anisotropic repulsion. 

As a first example we consider the model with = —1.96. Figure ^ (b) shows results for 
strip width L = 10 at different temperatures and constant chemical potential fi = —1.96. We 
have evaluated the coverage 6 = { rix^y ) = J2x,y f^x^y/M as well as the correlations 

Cd = ]^{nx,y{nx+i,yn+nx+i,y-i)) and Cx = {ux^yUx+i^y) . (2) 

These measure the probabilities of finding an occupied NN-pair {cx) or NNN-pair (c^) of 
Cd-atoms, i.e. the contribution of (2 x 1)- or c(2 x 2)- dominated regions in the system. 
Coverage and correlations can be obtained from proper derivatives of InZ^, or, as in the case 
of 9 and Cx , directly from the relevant eigenvector of the TM js], |^ . 

In addition. Figure |l| (b) displays results of Monte Carlo simulations of a system with 
M = 64 X 64 sites. In order to achieve reasonably fast equilibration we have applied a 
rejection- free algorithm [^], the results are in good agreement with the TM-calculation. In 
addition to the correlations (^) we determine order parameters which are associated with a 
perfect c(2 x 2) or (2 x 1) structure on one of the sublattices: 

-j^ xeven (x+j/) even 

x,y x,y 

Large values (< 9) of these quantities indicate long range order, whereas a homogeneously 
disordered occupation of the lattice would yield m2x2 = ^-2x1 = 6^/2. For the sake of 
breaking the sublattice symmetry, we have initialized the system with m2x2 = ^ for the 
equilibration dynamics. We have refrained from determining the order parameters within 
the TM-approach, which would require the introduction of additional staggered fields to the 
energy function (||). The TM-formalism offers a more suitable method to localize the phase 
transition [Q]. 

In the considered example, one observes a sudden drop of the coverage at T 0.3 when 
fi = —1.96 is held constant. Simultaneously, the system looses its long range order as indicated 
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Figure 1: 

a) Structural model of the c(2 x 2) and (2 x 1) reconstructions of the CdTe(OOl) surface 
1^, p. Shaded areas mark the corresponding primitive unit cells. Large filled circles represent 
Cd-atoms at the surface, open circles correspond to Te in the underlying half-layer, and 
small filled circles to the next, completed half-layer of Cd. Crosses represent empty sites in 
the simplifying lattice gas model. Note that the Te-atoms are dislocated according to the 
Cd-positions in the respective reconstruction. 

b) The phase transition at constant chemical potential ^ = —1.96 for Jx = —1.96. The lower 
panel displays results of the TM-calculation for L = 10 (solid lines) and MC~simulations 
(64 X 64 sites, single run): coverage 6 (triangles), correlations Cd (diamonds) and Cx (circles). 
The upper panel shows m2x2 (squares) and m2xi (crosses) for the same temperature range. 



by values m2x2 = '"T'2xi = in the simulations. This is also signaled in the properties of 
the relevant eigenvector in the TM-analysis The behavior is consistent with a first order 
transition, as it was investigated for similar models with isotropic or anisotropic interactions, 
see e.g. [|l|, ^, ^, ^, |5| and references therein. 

Here, however, also the NNN-correlation decreases rapidly at the coverage drop, while 
Cx displays a sudden increase and Cx > in the high temperature regime. This indicates that 
the phase transition also affects the short range correlations in the system: atoms order in 
rows of the (2 x l)-type without long range order. At = 1/2 the c(2 x 2) ordering is always 
preferred energetically. For significantly smaller coverages, however, the local rearrangement 
of atoms is possible and can be favorable if Jx ~ ^J^. Indeed, the degree of the prevalence of 
Cx over Cd depends strongly on the actual coverage as will be discussed below. 

We have followed the prescription outlined by Bartelt et al. |Q] for estimating the coverage 
discontinuity and phase boundaries for L ^ oo from three different strip widths. The results 
as obtained from L = 6,8, 10 are shown in Figure § for the models with Jx = —1.90 and 
Jx = —1.60, i.e. AE = 0.05 and 0.2, respectively. At low temperatures (III), an ordered 
phase with ~ 1/2 coexists with a disordered phase of low coverage. At higher temperatures, 
the system becomes homogeneously disordered (II) or ordered (I) depending on the coverage. 
For T — > oo, we expect the phase boundary (I/II) to approach the 9 = 1/2 axis. In this limit 
the infinite repulsion should be the only relevant interaction, columns of lattice sites decouple 
and the system is always disordered. This is in contrast to hard square models with isotropic 
NN-repulsion, where an extended regime (I) persists for arbitrary temperature. 

As an additional characteristics of the system we have determined the line T{9) where 
Cx = Cd and extrapolated for L — > oo. Right of the dashed lines in Figure ||, the c(2 x 2)- 
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Figure 2: Phase diagram of the model with Jx = —1.90 (left panel) and Jx = —1.60 (right 
panel), note the different temperature scales. Phase I is homogeneously ordered, in region II 
the system is homogeneously disordered, and in III the high and low coverage phases coexist. 
Solid lines represent the TM-extrapolation and symbols (circles) correspond to the results of 
Monte Carlo Simulations (M = 128 x 128) at constant coverage. The additional dashed lines 
(squares, respectively) indicate the values of (6',T), where Cx = Cd, hence it separates the 
region of c(2 x 2)-prevalence from the one where the (2 x 1) structure dominates. Statistical 
errors would be on the order of 0.05 for all the simulation data. 



structure is prevalent and vice versa. For small coverage, this characteristic line coincides 
with the boundary (II/III) of the coexistence region. Hence, for a range of coverages, the 
transition into disorder is accompanied by a simultaneous and discontinuous change of local 
ordering from c(2 x 2) to (2 x 1) arrangement of Cd-atoms. 

We obtain also a rough estimate of the phase diagram from additional Monte Carlo simula- 
tions at constant coverage. For this purpose, we apply a non-local algorithm which exchanges 
empty with occupied sites according to Kawasaki-like rates [ pC[| . The system is again initial- 
ized in an ordered c(2 x 2)-configuration for equilibration, and a rapid decrease of m2x2 with 
increasing T marks the transition into the homogeneously disordered phase. Figure ^ shows 
in both diagrams the results for M = 128 x 128, which are in good agreement with the 
TM-prediction. Within error bars, we obtain the same results by searching for a pronounced 
maximum in the fluctuations of order parameters, correlations, or energy. Note that this 
method is not suitable for detecting the transition into the homogeneously ordered region (I): 
simulations slow down considerably at almost maximal coverage and, furthermore, (I) and 
(III) become virtually indistinguishable in small systems. 

Figure |2| demonstrates the crucial role that the energy difference AE plays for the nature 
of the phase transition. With increasing AE, the tricritical point shifts to smaller coverage 
and higher temperature. Even more so does the line which separates c(2 x 2) from (2 x 1) 
prevalence. This feature might offer a qualitative explanation for the remarkable fact that 
the c(2 X 2)-(2 x 1) transition, which was investigated for CdTe in great detail, has not 
been found in ZnSe, so far. There, AE is expected to be significantly larger than for CdTe 
and the region of noticeable (2 x l)-dominance should indeed be smaller. Note that in the 
experimental investigation, integrated HRLEED-peak intensities provide information about 
local correlations, similar to Cx and Cd, rather than about long range ordering 

In summary, our model offers an interpretation of the c(2 x 2)-(2 x l)-transition in 
CdTe(OOl) as an equilibrium phase transition. At medium coverage the transition is, with 
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increasing T, from a coexistence regime into a homogeneously disordered phase. For smaU 
enough energy difference this phase transition is accompanied inevitably by a rearrange- 
ment of the vacancy structure from c(2 x 2)- to local (2 x l)-ordering. 

Of course, some of the detailed experimental findings cannot be accounted for in our simple 
model, see for instance Q for particular phenomena related to the relaxation of surface strain. 
For a more quantitative comparison with experiments, additional information is needed. A 
precise measurement of as a function of the temperature is difficult, but would reveal the 
path on which the system enters the (2 x l)-dominated region in the phase diagram. 

In a naive attempt to interpret our results quantitatively one would identify the dimen- 
sionless critical temperature (in units oi \ Jd \ = ^) with ~ 270°C, thus setting the scale 
for expressing the energy difference \2 + [/2 in physical units. For example, the model 
with Jx = —1.94 exhibits the desired transition with 9 ~ 0.4 at a temperature T ~ 0.3. This 
would translate into IS.E ~ O.OOSel^ which is significantly smaller than the value (0.03ey) 



given in |14, 15, 16]. DF-calculations yield /S.E at T = and the precise effect of higher 
temperatures on the relation of (free) energies is unknown. Furthermore, recent calculations 
have shown that the DF-results are very sensitive (up to a factor of about 2) to the number 
of atomic layers considered in the calculation ||2l|. Hence, a serious quantitative matching is 
not feasible unless more reliable estimates of Ai? become available. 

Another open question is, if and how our results for small values of 6 can be interpreted 
in the experimental context. Terminating layers of metal atoms with very low coverage are 
unstable in vacuum and the next (metal) layer is uncovered, see e.g. [0, 11 1. However, 
the presence of excess group VI atoms might stabilize an effective equilibrium situation with 
small metal coverage. As a test for this hypothesis we suggest to search for the structural 
transition of the ZnSe(OOl) surface under mildly Se-rich conditions. 

Our model also opens the possibility to study the shapes and sizes of domains, e.g. the 
regions of local (2 x l)-dominance in the disordered phase. Experimental data is available for 
the pronounced anisotropy of such domains Furthermore, we will study the equilibrium 
shape of isolated islands of atoms and its dependence on the temperature. This should allow 



for further comparison with experimental results as reported in [11|, for instance. 
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